Electron Transport in Disordered Graphene Nanoribbons 
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We report an electron transport study of lithographically fabricated graphene nanoribbons of 
various widths and lengths at different temperatures. At the charge neutrality point, a length- 
independent transport gap forms whose size is inversely proportional to the width. In this gap, 
electron is localized, and charge transport exhibits a transition between simple thermally activated 
behavior at higher temperatures and a variable range hopping at lower temperatures. By varying 
the geometric capacitance through the addition of top gates, we find that charging effects constitute 
a significant portion of the activation energy. 
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In recent years graphene has been celebrated for its 
potential as a new electronic material [H, Q. However, 
the absence of an energy band gap in graphene poses a 
challenge for conventional semiconductor device opera- 
tions. Previous work has shown that this hurdle 
can be overcome by patterning graphene into nanome- 
ter size ribbons or constrictions. The resulting trans- 
port gap formation can be most simply attributed to 
quasi 1-dimensional (ID) confinement of the carriers, 
which induces an energy gap in the single particle spec- 
trum jdl. Detailed experimental studies of disordered 
GNRs [3, H, i, E3, EJ 113, > however, suggest that this 
observed transport gap may not be a band gap. In an 
effort to explain these experimental results, various the- 
oretical explanations for the transport gap formation in 
disordered graphene nanostructures have been proposed, 
including models based on Coulomb blockade in a series 
of quantum dot s I13L A nderson localization due to edge 
disorder [13. [TBL \l(L 17l. HiL flit , and a percolation driven 
metal-insulator transition [20[ . Further systematic exper- 
iment is necessary to distinguish between these different 
scenarios. 

In this letter, we study the scaling of the transport gap 
in graphene nanoribbons (GNRs) with various widths 
and lengths at different temperatures. We find that four 
different energy scales can be extracted from transport 
measurement. From the scaling of these characteristic 
energies with GNR width and length, we find evidence 
of a transport mechanism in disordered GNRs based on 
hopping through localized states whose size is close to 
the width of the GNRs. 

GNRs with different lengths (L) and a widths (W) 
were fabricated following the procedures described in Q . 
Most experiments in this report were performed on back- 
gated GNRs on a substrate of highly doped silicon with 
a 285 nm thick Si02 gate dielectric. An example of 
such a device is shown in the inset to Fig. 1(a). We 
measured electron transport in a total of 41 GNRs with 
20 < W < 120 nm and 0.5 < L < 2/im at different tem- 



peratures 1.5 < T < 300 K. Additionally, we fabricated 
top-gated GNRs with 15 nm of hydrogen silsesquioxane 
(HSQ) and 10 nm of Hf02 as the gate dielectric material. 
The increased capacitive coupling allowed a comparative 
study of charging effects in back gated GNRs. 

GNR conductance is strongly suppressed for a region 
of back gate voltages Vg near the graphene charge neu- 
trality point 0, [111, IE IH [3, suggesting the for- 
mation of a transport gap. Fig. 1(a) shows low bias dif- 
ferential conductance G = dl/dV as a function of V g for 
a typical GNR. The transport gap region in back gate 
voltage, AVg, can be identified in this curve by extrap- 
olating the smoothed dG/dV g to zero [!, [Hj- We note 
that reproducible conductance peaks appear in the gap 
region 

IE BIH (left inset Fig. 1(a)), which are indicative 
of resonant conduction paths through localized states in- 
side the transport gap. In general, resonance peaks in 
the gap are less than 10 % of the G values outside of the 
gap region. 

The observed transport gap, AV g corresponds to an 
energy in the single particle energy spectrum: A m = 
hvF^/27rC g AV g /\e\, where vf = 10 6 m/sec is the Fermi 
velocity of graphene [23[ and C g is the capacitive coupling 
of the GNR to the back gate. This geometric capacitance 
is strongly dependent on ribbon dimensions and we cal- 
culate it using a finite element model, obtaining, for ex- 
ample, C g = 690 aF//im 2 and A m = 200 meV for the 
particular device in Fig 1. 

Away from the small resonant conductance peaks, the 
conductance is strongly suppressed in the transport gap, 
and the dominant charge transport can be described by 
thermally excited hopping between localized states [211 ] . 
We study the thermal activation of the off-resonant con- 
duction in this regime by measuring G m i ni the mini- 
mum conductance for a given sweep of gate voltage V g , 
at different temperatures (inset to Fig. 1(b)). Fig. 1(b) 
shows an Arrhenius plot for G m i n (T). Evidently, ther- 
mally excited transport exhibits two distinct behaviors 
at low and high temperature regimes, respectively, sepa- 
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FIG. 1: (a) Differential conductance of a GNR with W = 
36 nm and L = 500 nm, plotted as a function of back gate 
voltage. Dashed lines highlight measurement of AV g . Right 
inset shows an atomic force microscope image of the device. 
Scale bar is 500 nm Left inset shows a close-up of conduc- 
tance within the gap regime plotted as a function of V g — Vd, 
where Vd = 21 V is the gate voltage for the charge neutrality 
point, (b) T dependence of the minimum conductance of the 
same GNR in (a). The dashed line is a fit to simple activated 
behavior; the dotted line is a fit to variable range hopping 
with 7 = 1/2 and To = 460 K. An arrow highlights the po- 
sition of T*. 
temperatures. 



rated by a characteristic temperature T*. At high tem- 
peratures (T > T*), the transport is simply activated: 
Gmin ~ exp(— T^/2/c^T), where E a = 285 K is obtained 
from a linear fit of the Arrehenius plot (dashed line). 
At lower temperatures (T < T*), however, G m i n devi- 
ates from the simple activation behavior and decreases 
more slowly with decreasing temperature than the acti- 
vated transport would imply. In this low temperature 
regime, the overall behavior is consistent with variable 
range hopping (VRH), where G ~ exp(— (Xb/T) 7 ), with 
1/3 < 7 < 1/2 and a constant To, determined by the 
characteristics of the localized states [22| . 

The aforementioned GNR transport gap and temper- 
ature dependent characteristics are typical of all GNRs 
with W < 80 nm, so that A m , E a , and ksT* can be 
determined for each of these narrow GNRs. These three 
representative energy scales are plotted as a function of 
W in Fig 2. In this graph, we note that (i) there is a 
clear separation between these energy scales, setting a 
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FIG. 2: GNR transport energy scales : A m (solid), E a 
(shaded), and /c^T* (open) plotted as a function of GNR 
width. Circles correspond to ribbons of L — 500 nm. Trian- 
gles, squares, and stars correspond to ribbons of length 1, 1.5, 
and 2 /im, respectively. The dashed lines are the fits described 
in the text. 



general relation: A m > E a > /c^T* for given W; (ii) 
A m , E ai and T* depend sensitively on W but not T; 
and (iii) the energy scales are reasonably well described 
by inverse proportion to the lateral confinement of the 
GNR. The length independence can be noticed by com- 
paring characteristic energies of the GNRs with similar 
W but different L (represented by different symbols in 
Fig. 2), and suggests that these three energy scales are 
ID intensive properties of GNRs. To show this, we de- 
fine the normalized width w = (W — Wo)/ao? where 
ao = 0.142 nm is the carbon-carbon bond length and Wo 
is an offset introduced phenomenologically. Then, we find 
that all energy scales can be reasonably fit (dotted lines): 
A m = /SP m /w] E a = El/w] T* = T£/w with the propor- 
tionality parameters A^ = 36.3 eV, E® — 3.39 eV, and 
ksTo = 347 meV, respectively, with W = 12 nm held 
fixed for all three fits [241 ]. 

Edge disorder in the GNRs tends to induce wavefunc- 
tion localization, with a localization length that decreases 
rapidly with decreasing energy, resulting in a transport 
gap with strongly localized states at energies between the 
mobility edges [14| . The size of this mobility gap is larger 
than the clean band gap of an ideal ribbon; Querlioz et 
al. calculate the scaling prefactor A^ ^32.2 eV, aver- 



aged over many configurations of edge disorder 17]. The 
close match of our data to theoretical prediction sup- 
ports the view that atomic defects at the graphene edges 
create localized states. We point out, however, that the 
observed energy scales lie within the range of disorder po- 
tential fluctuation created by the charged impurities in 
the SiC>2 substrate [2(j, making it difficult to exclude the 
contribution of a substrate disorder induced transport 
gap, as discussed in a recent experiment on transport in 
thermally annealed GNRs [121 ]. 
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On the other hand, E^/A^ ^0.1; i.e., the activation 
energy at higher temperatures is an order of magnitude 
smaller than A m . This observation excludes the scenario 
that extended states carry current via thermal activation 
across the transport gap. Instead, we interpret the sim- 
ply activated behavior as a signature of ID nearest neigh- 
bor hopping (NNH) through localized states within the 
transport gap [19[. In this picture, disorder at the edges 
tends to produce a rapid variation in the local density of 
states over the whole width of the ribbon, blocking the 
conductive paths and leading to a quasi- ID arrangement 
of localized states [l6| . Martin and Blanter predict [l9[ 
that the energy spacing between nearest neighbor states 
is determined by ~ t'/w, where t' ~ 0.2t is the hopping 
matrix element between second nearest neighbor carbon 
atoms in graphene, so that E® ~ 2t' — 1.2 eV [19| . 
Our measured value for this scaling prefactor, 3.39 eV, 
is somewhat larger than this prediction, which may be 
explained by the contribution of a charging energy to the 
hopping energy E a , discussed in more detail below. 

The change of the transport behavior across the tem- 
perature T* allows a further comparison of our data to 
theory. In a very recent theoretical work, the NNH and 
VRH crossover is calculated to occur at T* = E a /k,BOL, 
where a ~ 8 was estimated numerically [25]. In our 
experiment, we obtain E^/ksT^ = 9.8, reasonably con- 
sistent with this theoretical prediction, lending further 
support to a model of charge transport via thermally ac- 
tivated hopping between localized states. 

An alternative approach to probing the GNR transport 
gap is measurement of the non-linear transport charac- 
teristics 0| . Fig. 3(a) shows differential conductance, 
dl/dVb as a function of V g and source-drain bias volt- 
age Vfe. Transport through the GNR at finite Vb shows a 
strong non-linear I — Vb characteristic when Ep is in the 
transport gap regime, which is most extreme when V g is 
near the charge neutrality point of the GNR (Fig. 3(b), 
black curve). The non-linear gap AVb can be defined 
where a steep increase of current appears in logarithmic 
scale (Fig. 3(b), green curve). 

In our previous study Q , the energy corresponding to 
eAVb was interpreted to be the band gap of the GNR. 
However, this naive interpretation should be carefully re- 
considered for edge disordered GNRs, where the charge 
transport is dominated by hopping through localized 
states. Indeed, from the plot of A Vb vs W (Fig 3(c)), 
we notice that AVb depends strongly on L, and is not 
well determined by W alone, unlike the previous three 
characteristic energy scales (A m , E a , and fe^T*). Since 
the charge transport in the disordered GNRs is diffusive, 
it is likely that electric field is driving transport in the 
transport gap. Indeed, if we convert AVb into the corre- 
sponding critical electric field £ cr = AVb/L, we restore a 
reasonable scaling behavior, where £ cr depends only on 
W and not on L (Fig. 3(d)). 

In disordered systems in which transport is domi- 
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FIG. 3: (a) Differential conductance as a function of V g and 
Vb measured in a GNR with L = 1 /im and W = 31 nm. (b) 
Current as a function of Vb with V g fixed in the off-resonant 
condition marked by the dotted line in (a). AVb is highlighted 
by the vertical dashed lines, (c) AVb as a function of W. 
Symbol shapes in (c) and (d) represent different GNR lengths 
following the convention set in Fig. 2 (d) The critical electric 
field Scr versus W converted from the data set in (c). 



nated by hopping through localized states, applied elec- 
tric field £ plays a similar role to temperature. Thus we 
can treat the electric field as an effective temperature: 
ksT e ff = eSL c , where L c is the averaging hoping length 
between localized states [27[. Noting that the transition 
from NNH dominated transport to VRH transport occurs 
at T*, we can estimate L c « ksT* / e£ cr . For most GNRs 
in this experiment we find that W < L c < 2W (Fig. 3(d) 
inset). The fact that L c > W supports our claim that 
hopping transport through the ribbons is effectively ID. 

Finally, we discuss the effect of Coulomb charging in 
GNRs. Several previous works have discussed the role 
of Coulomb blockade and charging effects on the trans- 
port gap in GNRs and graphene constrictions @, [H, [l3[ . 
In principle, in a GNR with hopping between localized 
states, we expect Coulomb interactions to open a soft 
Coulomb gap near the Fermi surface, which can be in- 
corporated into the total hopping energy E a in addition 
to the single particle energy level spacing t'/w, so that 
E a ~t'/w +_Ec, where E c is the Coulomb charging en- 
ergy. 



Til . 0, Q . In order to quantify the contribution of 
charging energy E c to the hopping energy E a , we perform 
a comparative transport measurement on GNRs with dif- 
ferent gate coupling. Fig. 4 shows the temperature de- 
pendent minimum conductance G m i n (T) for a back gate 
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FIG. 4: Temperature dependence of the conductance mini- 
mum for dual gated (circles) and back gated (triangles) GNRs 
with the similar W and L. The dashed lines are Arrhenius fits 
in the high temperature regime. The inset shows SEM images 
of back gated (left) and dual gated (right) devices. Scale bar 
represents 500 nm. 



only GNR (device I) and a GNR with both top and back 
gates (device II) with the similar W and L. While device 
I has usual capacitive coupling to the back gate, (i.e., 
C 1 w C g ), C 11 for device II is much closer to the top 
gate, leading to a larger capacitance: C 11 /C 1 « 4. From 
the thermally activated Arrhenius behaviors in the high 
temperature regime (dashed lines), we obtain the acti- 
vation energies of the two devices, E T a — 15 meV and 
E^ 1 = 8.4 meV averaged over two devices of type I and 
four of type II. Considering the smaller charging energy 
contribution for a top gated device, smaller values of the 
activation energy are indeed expected, if Coulomb effects 
are appreciable in the GNR. 

Employing the ratio E^/E^ 1 « 0.5, we now can es- 
timate the charging energy contribution quantitatively. 
Assuming that the single particle energy level spacing 
t'/w is similar for both GNRs due to their similar di- 
mensions, we obtain E* 1 — E]} = El — E]. = t'/w, where 
the charging energy ratio of device I and II are given 
by El /El 1 = C II /C I « 4. The resulting estimate for 
the charging energy contribution, El /El w 0.6, indicates 
that the Coulomb charging effect provides a substantial 
portion of the activation energy. 

In conclusion, we investigate length and width depen- 
dent resistance scaling in GNRs. Temperature dependent 
and electric field dependent transport characteristics in- 
dicate that charge transport in the transport gap of the 
disordered GNR is dominated by localized states, where 
the Coulomb charging effects play an important role. 
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